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Abstract 

This M.Sc. thesis in Engineering Physics is an overview of the 
present theory of quark-gluon plasma (QGP) as well as an analysis of 
the stability criterion for possible stable cosmic QGP objects left over 
from the quark-hadron transition in the early Universe. It covers fun- 
damental ideas of the formation and decay of the plasma, including the 
standard model, QCD, and the MIT bag model. I discuss the equation 
of state of a QGP and the possible signatures for a plasma created in 
heavy-ion collisions. Astrophysical aspects of QGP are put forward, 
including compact stars and the quark-hadron transition in the early 
Universe. The possible role of QGP objects as cosmic dark matter 
is mentioned. The analytic part is an investigation of possible stabil- 
ity among cosmic QGP objects from the early Universe. A model is 
suggested where a pressure balance makes a QGP stable against gravi- 
tational contraction and hadronization. The mass/radius relationship 
for stability also forbids a direct gravitational collapse. Finally, the 
possibility of stable cosmic QGP objects is critically discussed. 
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Introduction 



This thesis is devoted to the subject of quark-gluon plasma, QGP, a fairly 
new branch of high-energy physics. A lot of effort has been put into the 
area, primarily into experimental detection of QGP production in heavy-ion 
collisions, but as of today, no undisputed signal seen. 

The largest experiments up to date in the search for a QGP is at the SPS 
accelerator at CERN in Geneva, where physicists have been looking for a 
QGP for several years without success. These experiments are still running, 
and will hopefully be continued in one form or the other when the new large 
hadron collider, (LHC), currently under construction at CERN, will be taken 
into use around the year 2005. 

In the US, a major effort to produce and detect a QGP is underway at the 
Brookhaven National Laboratory, where a new relativistic heavy-ion collider, 
(RHIC) is being built for a planned operation by 1999. 

This thesis is not going into details about the experimental efforts to 
detect a QGP. I will instead focus on the theory of a QGP together with 
an excursion into the exciting field of astrophysics. The thesis is split up 
into two parts. One is an overview of the physics of a QGP, including a 
short introduction to quantum chromodynamics (QCD) and the so-called 
MIT bag model, and the other describes the results I have achieved about 
some astrophysical aspects of QGP. 

In the overview I have tried to include those aspects of QGP that are 
widely stated as fundamental and well established. It starts with a short 
overview of the standard model and describes the basic properties of quarks 
and gluons. The gauge theory (QCD) describing the strong interaction is 
given some attention as well as the QCD-inspired MIT bag model of quarks 
confined in hadrons. Since the confinement mechanism is crucial in the phase 
transition from hadronic matter to QGP, a rather extensive description is 
given, as well as a discussion of the phase transition itself. 

The equation of state of a QGP is given, i.e., the relation between the 
pressure and temperature inside the plasma. This is followed by a discussion 
of the formation and decay of a QGP, which leads into a subsequent overview 
of possible experimental QGP signals. 
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Various astrophysical aspects of QGP are discussed, and the role of QGP 
in the phase transitions in the early Universe is mentioned, as well as the 
current theory governing compact stars, i.e., neutron, quark and hybrid stars. 
This provides a link between the QGP overview and my own calculations and 
results. 

I examine the possibility of stable cosmic QGP objects surviving from 
the quark-hadron transition in the early Universe. I show, within the theory 
of general relativity and the equation of state given by QCD, that stable 
configurations can occur when the size and the mass of the QGP have a 
certain relationship derived from the so-called Tolman-Oppenheimer-Volkoff 
equation. 

At the end I discuss some general aspects of QGP and point out certain 
discrepancies between the model emerging from heavy-ion collisions and ideas 
applied to astrophysical QGPs. 

In two appendices, chiral symmetry and the physics of phase transitions 
are discussed. 
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Chapter 1 
Quarks and gluons 



1.1 The standard model 

Most particle physicists today believe that the standard model of elementary 
particle physics more or less describes the fundamental building blocks of 
matter along with their interactions (apart from gravity). It comprises the 
theories of the electroweak and the strong interactions. The standard model 
tells us that we have two groups of elementary particles, leptons and quarks. 
In addition, the different types of interactions included in the model are due 
to exchange particles, in the form of vector bosons. 




The leptons can only interact by the electroweak interaction, i.e., the uni- 
fied electromagnetic and weak interaction. They do not feel the strong force 
mediated by the gluons. The quarks, on the other hand, interact strongly, 
weakly and electromagnetically. 
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Leptons and quarks obey certain empirical particle-number conservation 
laws. If the neutrinos are massless, one could speak of conservation of lepton 
type, i.e., conserved electron (v e ,e), muon and tau {u t ,t) numbers. 
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All leptons have antiparticles with opposite electric charges and lepton 
numbers. 

Quark flavour is conserved by the strong interaction, but not by the weak 
interaction. One can speak of a total quark number due to the stability 
of the proton. Up to this date, experiments and observations are consistent 
with total conservation of overall quark number and of lepton type, but these 
conservation laws are not a consequence of any known dynamical principle. 
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All quarks have their antiparticles, the antiquarks. They have the same 
quark quantum numbers and electromagnetical charges apart from a minus 
sign. 

The three types of interactions included in the model are all mediated by 
exchange of vector bosons. The mediator of the electromagnetic force is the 
photon, those of the weak force are the W ± and the Z°, and the strong force 
is mediated by the gluons. The photon carries no electric charge and there 
is no direct interaction between photons. The gluons actually carry colour 
charge, a sort of "strong" charge, and therefore interact with each other. The 
gluons and the photons are presumably massless, but the W ± and the Z° are 
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quite heavy. Therefore, the range of the weak interaction is very short, about 
1CT 18 m. The strong interaction has a limited effective range of 10~ 15 m due 
to colour-screening effects and gluon self-interaction. 



Particle 


Mass (GeV/c) 
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z° 


91.187 ±0.007 


7 


< 6 x 10~ 25 



1.2 Quarks 

Experimentally, one has found evidence for six quarks; up (w), down (d), 
charm (c), strange (s), top (t) and bottom (6), which all are fermions. Their 
masses differ much, from the up with a mass of around 300 MeV/c 2 to the 
top-quark with a mass around 180 GeV/c 2 [1]. All quarks have electric 
charges; ±|e for the u,c and t and — |e for the d,s and b. Each quark is said 
to represent a separate flavour. 

Quarks, as fermions, obey the Pauli principle, which presented some ma- 
jor difficulties when the A ++ resonance was discovered. The A ++ resonance 
is a spin/parity J p = | + particle consisting of three up quarks with parallel 
spins: 

|A ++ >= |wVw T > . 

This situation is not favourable in any way. Three identical particles 
with parallel spins violate the Pauli principle. This problem was solved by 
introducing a colour charge [2]. The colour quantum number can take three 
basic values, e.g. "red", "green" and "blue". Anti-particles can be "antired", 
"antigreen" and "antiblue" . The Pauli principle is now obeyed, provided that 
the three quarks in a baryon have different colours. The A ++ resonance has 
an antisymmetric wave function: 

|A ++ >= \u\u\ul > . 

When the quark concept was invented in the mid 1960s the physicists were 
able to categorize and describe most that had been discovered in accelerator 
experiments. The particles that are made of quarks were called hadrons; 
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strongly interacting particles. Hadrons with half-integer spin were called 
baryons and hadrons with integer spin mesons. 

The proton and the neutron are two examples of hadrons. Each consists 
of three quarks, 



and have electric charge +e and 0, respectively. In the proton/neutron case, 
the three quarks that give the nucleon its properties are called valence quarks. 
Virtual quark-antiquark pairs are continuously created and annihilated inside 
the nucleon. 

All in all, the quark family consists of six quarks, with different flavours, 
and they all carry electric charge as well as colour charge. 

No free quark has ever been detected, in spite of several extended searches. 
Therefore one believes that quarks are confined inside hadrons, and that the 
strong force potential between the quarks increases as the distance gets larger. 
One would then need an infinite amount of energy to separate two quarks 
from each other. Hence hadrons are colour neutral objects. The quarks 
forming a hadron must have a colour combination rendering a colourless, 
"white" particle. For instance, the spin-0 meson tc + has three possible colour 
configurations. 



A physical n + is equal to a quantum-mechanical mixture of these states. 



The force that binds quarks together is the strong interaction. It is mediated 
by its exchange particles, the gluons. The gluon is a spin-1 particle with no 
mass, and hence the strong interaction would have a long range would it not 
be screened. The effective range of the force is of the order of 10 -15 m. The 



\p >= \uud > 



n >= \udd >, 




1.3 Gluons 
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Figure 1.1. Graphical picture of the relation between the different colours and 
two examples of colourless, white states. 

gluons carry both colour and anticolour. Since 3x3 colour combinations 
exist, the gluons form two multiplets of states, an octet and a singlet. It 
is possible to construct all colour states from the octet and therefore there 
are eight different gluons. The ninth state is the totally symmetric state 
RR + BB + GG, which is colourless, and therefore plays no role in the strong 
interaction. 

Since gluons themselves carry colour they can interact with each other. 
This is the main difference between the strong and the electromagnetic in- 
teraction, since photons carry no electric charge. Gluons can be emitted and 
absorbed, created and annihilated, when a strong interaction is involved in 
a process. 

Since quarks make up the proton and the neutron, gluons mediates the 
force that keep the proton and the neutron together. The strong nuclear 
force that makes the nucleons in a nucleus stick together, is in quark-gluon 
terms a leakage of the forces that keep the nucleons together, e.g. in the form 
of quark-antiquark pairs. 

Since all quarks have the same possibility to carry a certain colour regard- 
less of flavour, the strong interaction is flavour independent, i.e., all sorts of 
quarks have identical strong interactions. 
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Figure 1.2. The fundamental interaction Feynman diagrams of the strong 
interaction. From left to right, splitting of a gluon into a quark- antiquark 
pair, emission of a gluon, and two gluon self- couplings. 

One consequence of the properties of the strong interaction is that the 
strength of the strong force decreases when the quarks are close to each other. 
This property gives the quarks what is called asymptotic freedom at small 
inter-quark distances. 

The strength of the strong interaction is described by the strong "run- 
ning" coupling constant a s . It has a direct dependence on the squared four- 
momentum transfer in the quark process: 



N/ is the number of quark flavours involved and A is a scaling parameter, 
A = 0.2 ±0.1 GeV. From this expression, it is easy to see how a s decreases 
with increasing momentum transfer. This is equivalent to decreasing a s when 
the distance to a quark decreases. Hence, asymptotic freedom and colour 
confinement are implied by the expression for the strength of the interaction. 

A more detailed account of the quantum field theory of the strong inter- 
action (QCD) is given in the next chapter. 




127T 



(1.1) 



(33 - 2N f )ln(Q 2 /A 2 ) ' 
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Chapter 2 

Quantum chromodynamics 



2.1 The concept of colour 

A new quantum-mechanical gauge theory of quark interaction was born in 
the 1960s and 1970s. It was called QCD, in analogy with quantum electro- 
dynamics (QED). 

The basic idea with QCD is invariance against arbitrary rotations in 
colour space. Since the complex rotations of an arbitrary vector in three- 
dimensional space are described by unitary 3x3 matrices of unit determinant, 
the symmetry group of the gauge transformation is SU(3). 

When the concept of colour was invented, the main reason was to make 
the quarks obey the Pauli principle, but at that time there was no evidence 
that there are exactly three colours. The experimental verification of the 
number of colours came in the beginning of the 1980s [3]. The ratio of the 
production of hadrons from e + e~ collisions to that of muon pairs reveals 
an interesting fact. This is based on an assumption that the production of 
hadrons proceeds through creation of a quark- ant iquark pair, which subse- 
quently fragments into hadrons. In that case, the ratio of the cross-sections 
can be related to the sum over the square of the electromagnetical charges of 
all sorts of quarks that can be created. Defining N c as the number of colours, 
one gets 
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for i = w, c/, s, c, 6 quark types. The experiments yield N c = 3 without doubt 



The principle behind the QCD Lagrangian is the demand for invariance un- 
der local colour rotations, i.e., the gauge matrix U changes from one point in 
space to another. In quantum electrodynamics (QED) the governing princi- 
ple of gauge transformation is changes in the phase of the wave function, as 
described by the one-dimensional U{\) symmetry group. 

The approach in deriving the QCD Lagrangian is the same as in QED 
modified to fit the gauge invariance of the strong interaction. The funda- 
mental QED Lagrangian is (see, e.g. [4]): 



where F^ v = d^A u — d v A tl is the electromagnetic field strength tensor, and 
An is the electromagnetic vector field. 

Due to the colour property, the wave function of a quark has three com- 
ponents in colour space, ip = ijp ri ip gi ipb)- A colour gauge transformation, 
ip — > Utp, is described by a unitary 3x3 matrix U with det(U) = 1. U can be 
written as the imaginary exponential of a Hermitian matrix L, U = exp(iL), 
where L* = L and tr(L) = 0. All traceless Hermitian 3x3 matrices can be 
expressed as linear combinations of the eight A-matrices [5] 



[3]. 



2.2 The QCD Lagrangian 



Lqed = i^id^ + ieA^ - miH) - -F^F, 



(2.2) 




(2.3) 
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[A a , A 6 ] = 2if abc \ c and [A a ,A 6 ] 



+ 



= -5 ab + 2d abc \c- (2.4) 
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These so-called Gell-Mann matrices, or |A a to be more exact, are the eight 
generators of the Lie group 577(3), and f a bc and d a bc are the antisymmetric 
and symmetric structure constants. To make the rotation U space-dependent 
the real parameter 6 a must vary in space, 6 a = 9 a {x). This leads to 



U (x) = exp 
Thus ip — > U(x)ip gives 



1 8 



2 

a=l 



(2.5) 



d^Uix)^] = U8^ + (^C/)V = U[d^ + U*{dJJ)il>\. (2.6) 

Since Ud^ip should be present in the right-hand side one has to introduce a 
colour potential, the analogue to the electromagnetic potential in QED. 
That potential is called A^; a 3 x 3 Hermitian matrix. This A^ can be 
represented as a linear combination of Gell-Mann matrices: 

1 8 

A,{x) = -Y,Al{x)\ a . (2.7) 

Z o=l 

The field A^(x) is called a Yang-Mills field [6], and if A^(x) changes during 
a colour rotation according to 

A^U^U-i-U^U), (2.8) 

the derivative (d^ — igA^) will be invariant under such gauge transformations. 
Here g is the so-called strength factor of the strong field. Since the basic 
relation is the Dirac equation one can choose a potential that leaves the 
derivative invariant under such a gauge transformation making the whole 
equation invariant. The replacement defines a covariant derivative, = 

The next step is to add a kinetic term to the Lagrangian, and a first 
guess is a QED-like term, —jF^F^. It turns out that one has to modify 
the definition of the field strength tensor F^ in order to keep the theory 
gauge invariant. Since no interactions take place between two photons, QED 
is an abelian theory and the generators of the U(l) group commute. QCD, 
on the other hand, is a non-abelian SU (3) local gauge theory. The generators 
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are non-commutative, due to the fact that interactions take place between 
the gluons. The field strength tensor therefore changes to 

F; u = d,A a v - d v Al + gf^Al, (2.9) 

and this definition makes it form invariant under a local colour gauge trans- 
formation, F^y — > U*F^ U U. 

The complete QCD Lagrangian is then 

Lqcd = - igA^ - - \F a ^ . (2.10) 

Notice the resemblance with the QED Lagrangian, eq. (2.2). All the dif- 
ference lies in the self-interaction term in the definition of F^ and in the 
exchange of +ieA^ — > —igA^. 
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Chapter 3 



The confinement of coloured 
particles 

3.1 The confinement mechanism 

The non-observation of single quarks led the physicists to postulate that no 
particle can exist in a coloured state, and as a consequence, to establish the 
concept of quark confinement. 

QCD has complex field equations that make it unsuitable for "exact" cal- 
culations. This leads to extensive use of perturbative methods. These are, 
on the other hand, valid only in certain regimes of the perturbative expan- 
sion parameter. In QCD, this parameter is the strong coupling constant, 
a s . It gets small when the four- momentum transfer is large or, equivalently, 
the distance involved is small as illustrated by eq. (1.1). That makes per- 
turbative QCD valid only in processes that have these characteristics. This 
is unfortunate because the distance when a s gets too large for a perturba- 
tive expansion, is about the same as the dimensions of a hadron. Therefore, 
the confinement mechanism cannot be derived directly from the QCD La- 
grangian, and confinement cannot even be proven. Instead, one has to use 
"QCD-inspired", phenomenological models to describe confinement. 
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3.2 The QCD vacuum 



The expression for the strong coupling constant eq. (1.1) is derived from the 
full gluon propagator, D(q 2 ) where q 2 = —Q 2 . This propagator is essentially 
the boson propagator, but it has been modified due to vacuum polarization. 

It is possible, to reach the same quantative conclusion involving the QCD 
vacuum, which follows a line of reasoning by Gottfried and Weisskopf [7]. 

Let us consider the empty vacuum, i.e., one where there are no gluons. 
To this vacuum, add a pair of gluons of opposite colour charges and spin- 
components, and with an average separation r. The energy for this pair 
is 



where A and C are constants. The first term is the kinetic energy and the 
second one is the potential energy. r is the distance where a s gets "large". 
For r <C ro, a s is small and E(r) positive, but near ro, E(r) < 0! When 
r > ro the expression is no longer valid, but it should suffice to assume that 
E(r) increases with r. 




for 



(3.1) 



E 




r 



Figure 3.1. The energy E of the gluon pair versus distance r. 
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If this picture is correct, it seems to be energetically favourable to create 
a gluon pair with opposite spins and colours out of the "empty" vacuum, to 
take advantage of the negative energy when r ~ r . 

The true vacuum can now be described as follows. The "empty" vacuum 
is unstable. There exists a state of lower energy that consists of cells, each 
containing a gluon pair in a colour- and spin-singlet state. The size of these 
cells is of order r . 

This give one a crude description of how a colour-neutral assembly of 
quarks is immersed into the gluon vacuum. The gluonic cells will be displaced 
by the quarks, and therefore the quarks will find themselves in a gluon-free 
"bubble" or "bag". A state with quarks and gluons in the same bag has a 
higher energy, and therefore the quarks "push away" the gluon field. The 
size of this "bag" is of the order of 1 fm. Since the true vacuum has an energy 
density lower than the "empty" vacuum, the energy density in the bag must 
be positive. Actually, it is proportional to the bag volume, V: 

E B = BV. (3.2) 

B is called the bag constant and has the dimension of (energy) 4 provided 
that one sets c = h — 1 so that distances are counted in 1/eV. The region 
outside the bag exerts a pressure on the bag which is counteracted by the 
kinetic energy of the quarks. It can be shown [7] that this picture describes 
the observed structure, size and low-lying spectra of hadrons reasonably well 
if 5 1 / 4 m 150 MeV. The pressure on the bag then amounts to ~ 10 23 atm. 

Now one can compare the field of two opposite electric charges with the 
field of two opposite colour charges at a separation larger than r . In the 
electromagnetic case the field lines spread when the separation grows. The 
number of electric field lines crossing a unit area decreases like r~ 2 . In the 
colour case, the pressure of the true vacuum compresses the field lines into 
a tube of diameter r®. When r ^> tq the number of field lines per unit area 
within this tube remains constant, leading to a constant force, i.e., a linear 
potential: 

4>{r) = ar for r r . (3.3) 
The experimental value of the constant a is a ~ 0.6 GeV/fm. 
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Figure 3.2. Phenomelogical picture comparing the nature of the strong force 
with the electromagnetical force. 



In conclusion, there is a constant force acting on a quark if one tries to 
remove it from a hadron to distances larger than 7*0. If one would like to 
remove it completely one would need an infinite amount of energy. However, 
long before that, the colour field would "break up" into new quark- ant iquark 
pairs. 



Since the potential in QCD grows indefinitely with distance, coloured parti- 
cles are confined to each other. The so-called bag models have been invented 
to take this into account, because the exact QCD formalism is virtually im- 
possible to solve for bound hadronic states (just like QED for atoms). In 
these models, quarks are confined to a certain hadronic volume V. The 
Dirac equation for the quarks becomes 



Then, the quark current through the surface is zero. This makes the bound- 
ary condition read 



3.3 Bag models 



i^d^ip — mip = inside 
ip = outside. 



(3.4) 



(3.5) 
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This bag model was proposed by Bogolioubov in 1967 [8]. 

3.4 The MIT bag model 

In the mid 1970s a group of physicists at the Massachusetts Institute of 
Technology (MIT) showed that Bogolioubov's bag model leads to energy- 
momentum conservation violation at the bag surface, unless the internal 
pressure from the interior of the bag is balanced by an external pressure. This 
led to a modified model, the MIT bag model [9-11], where the true vacuum 
inside the bag is partially destroyed by quarks carrying colour. This mixture 
of quarks and true vacuum made it possible to treat the physics of the interior 
of the bag by perturbative QCD, a so-called perturbative vacuum inside the 
bag. This change in the bag model leads to a new boundary condition. The 
requirement of pressure balance at the surface is written as 

i 

where B is the bag constant, and the sum runs over all quarks contained in 
the bag. For a spherical bag this condition is equivalent to the requirement 
that the total energy contained in the bag volume should be at a minimum 
with respect to the bag radius R. Hence ^ = 0, where 

M(R) = + ^-BR 3 , (3.7) 

Xi = u>iR, and u>i are the eigen-frequencies of the solutions to the Dirac 
equation. The equilibrium radius is obviously 

It is possible to include real gluons in the bag even without quarks. Such 
states, "glueballs" have not been found, although signals of hybrid glue- 
ball/quark states are sometimes reported. The appropriate boundary con- 
dition for the colour field is obtained from the requirement that the colour- 
electric field should not be able to penetrate outside the bag. The situation 
is analogous to the boundary condition in classical electrodynamics for a 
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medium with e = and [i — oo. By virtue of Gauss' theorem an integra- 
tion of the colour-electric field yields the result that the total colour charge 
contained in the bag volume vanishes. 

There are several flaws in the MIT bag model, for instance, 

• that the mass of the pion comes out too large in all versions of the 
model 

• that the boundary condition is not chirally invariant, or equivalently, 
the MIT bag model violates PCAC, Partially Conserved Axial Vector 
Current. 



Several attempts have been made to overcome these difficulties. In the 
chiral bag models a pion field has been added as an independent degree of 
freedom [12,13]. In the cloudy-bag model [13,14], or in the Tel Aviv model 
[15], the pion cloud is allowed to penetrate into the bag. In these models, 
the main goal has been to make the model chirally invariant at the surface. 

One of the models, the soliton bag model, where the large mass of a 
quark outside the bag is generated by the coupling to a scalar field, is very 
well suited for the description of the dynamical properties of the bag, such 
as vibrations of the bag surface. This model could be useful in the study of 
giant QGP bags. 

It is important to remember that all the bag models are only attempts to 
construct a model that explains the behaviour of the particles as measured 
in experiments. As of this date, the mechanism of confinement is well es- 
tablished experimentally, but since perturbative QCD cannot be used in this 
region the theoretical understanding of confinement is incomplete. 
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Chapter 4 



The transition to a quark-gluon 
plasma 

4.1 Hadrons as systems of quarks 

Ordinary matter is made out of leptons and hadrons. Hadrons are strongly 
interacting particles consisting of quarks. We know from statistical physics 
and thermodynamics that macroscopic systems of ordinary matter exhibit 
changes in phases when the environment changes in a specific way. Why 
should not nucleons, built of quarks, exhibit some change in phase when the 
environment changes? 

It is experimentally justified to describe collisions of nuclei, i.e., systems 
of hadrons, in thermodynamic terms. The predicted temperatures in the 
collisions are in accord with experimental results. This temperature can be 
estimated from the energy distribution of the emitted fragments. The only 
discrepancy appears when the impact energy gets very high, the temperature 
does not increase as fast as the model suggests. The temperature seems 
to approach some kind of plateau when the collision energy grows. The 
temperature of boiling water does not change during the phase transition 
when the water enters its vapour phase. Would colliding nuclei approach 
some kind of analogue phase change? It turns out that the levelling-off of 
the temperature is due to the fact that higher-mass hadrons are created in 
the collision. The energy in the collision is high enough to convert into mass 
and create heavier hadrons than protons and neutrons. 
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It is possible to look at such a creation of hadrons as a kind of phase 
transition. It has a latent heat just as the liquid-vapour transition in water. 
The temperature in the collision remains the same until the most massive 
hadron have been created. Then the temperature starts to increase again. 
Hagedorn has suggested that there exists such a limit, and the limiting tem- 
perature should be about 1.5- 10 12 K [16]. This is not so far away from current 
accelerator experiments. 

4.2 From hadronic matter to QGP 

The ultimate phase transition would be the transition from hadronic matter 
to a quark-gluon plasma. In the plasma phase the nucleons and the higher 
mass hadrons created in the collision lose their identity as individual particles. 

It is important to remember the difference between the bonding of an 
electron to a nucleus and the bonding between quarks in a nucleon. It is 
possible with a finite investment of energy to separate an electron from a 
nucleus, and turn a gas into a plasma where all electrons move freely. But in 
a nucleon, one would need an infinite amount of energy to separate a quark 
from its environment. 

There are two ways to create a plasma in a gas of atoms. One is to heat 
the gas until the collisions between the atoms are violent enough to rip away 
all electrons. Another is to compress the gas, driving the atoms into close 
contact and make their electronic shells overlap. Under these conditions the 
electrons are deconfined and can move freely from one atom to the other. 

As explained earlier, the quark interaction is negligible at short distances 
(asymptotic freedom). On the other hand, when the distance increases the 
interaction also increases (confinement). These characteristics rule out the 
first way (heating) to create a quark-gluon plasma, but not the second one 
(compression). Hence, when making a quark-gluon plasma one does not need 
to set the quarks free, only to push them into the same bag. It should be 
added that it is at least theoretically possible to create a plasma by heating, 
although the phase transition is not obtained by quark deconfinement but 
by a more subtle mechanism [17]. 
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4.3 The QCD phase transition 



At low energies, the QCD vacuum is characterized by nonvanishing expecta- 
tion values of vacuum condensates, the quark condensate having < ipip >^ 
(235 MeV) 3 and the gluon condensate having < a 8 G llv G' J,v >^ (500 MeV) 4 . 
The quark condensate describes the density of quark- ant iquark pairs found 
in the QCD vacuum and is the expression of chiral symmetry breaking 1 . The 
gluon condensate measures the density of gluon pairs in the QCD vacuum 
and is a manifestation of the breaking of scale invariance of QCD by quantum 
effects. 

Broken symmetries in nature are often restored through phase trans- 
formations at high temperatures. Some examples are ferromagnetism, su- 
perconductivity and the transformation from solid to liquid. In the quark- 
gluon plasma case, the phase transition from hadronic matter to quark-gluon 
plasma would restore the broken chiral symmetry. Nuclear matter at very 
high temperature would exhibit neither confinement nor chiral symmetry 
breaking. 

It is known from quantum field theory that chiral symmetry is broken 
if the particles involved have non-zero mass. But in the QCD vacuum, a 
non-zero value of the quark- ant iquark condensate expectation value < ipip > 
has the same effect. The density of spin-0 quark- ant iquark paris has the 
chiral decomposition qq = q L qn + qi/lR, where subscripts R ("right-handed") 
and L ("left-handed") denote states with spin parallel and antiparallel to 
the direction of motion. Hence the broken vacuum state contains pairs of 
quarks of opposite chirality. If a left-handed quark annihilates a left-handed 
antiquark, then the process is perceived as a change of chirality of the free 
quark, exactly the same effect as a non-vanishing quark mass. In reality 
though, the light quarks have mass and the chirality of a light quark is not 
conserved, even if the quark condensate vanishes. 

When the temperature increases the interactions occur at shorter dis- 
tances, mainly by the weak coupling because the long range strong interaction 
gets dynamically screened. Finite temperature perturbation theory shows 
that the strong coupling constant a s falls logarithmically with increasing 

1 See Appendix A for an explanation of chirality and chiral symmetry. 
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temperature [18,19]. The result that the quark-condensate parameter van- 
ishes at high temperature, makes it likely that the transition between the 
low-temperature and the high-temperature manifestations of QCD shows a 
discontinuity, i.e., a phase transition. The order of this transition is believed 
to depend on the number of light, quark flavours. Two massless flavours 
would generate a second-order transition [20,21] while three massless flavours 
give a first-order transition. 

Numerical simulations indicate that the transition temperature lies in the 
range 150 ± 10 MeV [22]. 
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Chapter 5 

The properties of a QGP 



5.1 Equation of state 



Considering the energy scale in the hot QGP, i.e., ~ 200 MeV, the approx- 
imation that the plasma only contains u and d quarks which are massless 
should be quite reasonable. With this approximation and the assumption 
that one can neglect all quark interactions inside the plasma, one can derive 
the equation of state [4]. 

First, one needs to count the degrees of freedom for the constituents: 



The energy density residing in each degree of freedom is calculated separately 
for the quarks and the gluons. The gluons form, without the interactions, an 
ideal relativistic Bose gas of temperature T which gives the energy density 



where (3 = 1/T. For the quarks and the antiquarks one has to introduce 
a chemical potential /i because, in general, there will be a slight surplus of 
quarks over antiquarks in the QGP created from normal atomic nuclei (heavy 
ions). At zero temperature, \i is the energy needed to add another quark to 
the plasma. Since no antiquarks are present at T = the energy necessary 
to add an antiquark is zero. This does not imply that ji = 0, because 



Gluons: N g = 2(spin) x S(colour) = 16 



(5.1) 



Quarks: N q = 2(spin) x 2( flavour) x 3(colour) = 12. 



(5.2) 




(5.3) 
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the additional antiquark may annihilate a quark and release the energy /i, 
assuming that the quark lies at the surface of the Fermi sea. The chemical 
potential of the antiquarks must therefore be chosen as —fx. 

The energy density for the quarks, treated as a relativistic Fermi gas, 
cannot be calculated analytically if /x, T ^ 0. However, the energy density 
for a quark and an antiquark is 

r d 3 p p 

q " J (27r) 3 [e^) + l]' 1 ' 

/d^p p 
(2tt) 3 + 1] ' (5 ' 5) 

/d 3 p ( p P i 

J^Y { + i] + pra + 1] J ^- 6 > 



120 4 8tt 2 

Considering baryon-number symmetric quark-gluon matter, i.e., fx — 0, and 
multiplying with the respective degrees of freedom, the energy density be- 
comes 

E = 16£ 9 + 12( £ , + E.) = ^T* = ( M ^ v )* [GeV/W] . (5.7) 

The expressions above use = c = h = I, giving T the dimension of energy, 
and length the dimension of (energy -1 ). The energy density inside a nucleon 
is four times the MIT bag constant, E N = AB m 300 - 500 MeV/fm 3 , which 
together with a transition temperature of ~ 150 MeV gives an energy density 
of the QGP of ~ 1 — 1.5 GeV/fm 3 . From this one can see that one needs at 
least a factor of two higher energy densities than inside a nucleon in order to 
make a QGP. 

But what if fx ^ 0? One can still use eq. (5.7) but fx must now be 
computed, i.e., one has to know its relation to the baryonic density This 
density is one third of the density difference between the quarks and the 
antiquarks multiplied with the degrees of freedom: 

/dp 
eflP-zO + 1 ' (5 ' 8) 
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dp 



(5.9) 



n„ 



n,= -T 



67T 2 



n B = A(n q 



n-) = —T 2 + 



(5.10) 



Hence n B = where E = E q + E q is taken from eq. (5.6). The pressure 
and entropy of the plasma are given by 



p 



3 ' 
IdE 
3df' 



(5.11) 
(5.12) 



Under what conditions should the bag containing the plasma be stable? A 
fair assumption would be that the external vacuum pressure, characterized by 
the bag constant B, should not exceed the internal pressure, i.e., p = j > B. 
The critical temperature T c and the critical chemical potential /i c can be 
calculated using the eqs. (5.7) and (5.11) evaluated at the critical pressure 
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Considering also perturbative interactions in the QGP, a second-order 
thermal perturbative calculation gives [4] 
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where / is the quark flavour (neglecting quark masses). This formula includes 
also the effect of strange quarks (ss pairs), provided they can be treated as 
massless. However, it does not include short-range spin effects in QCD, so- 
called colour-magnetic forces. These tend to favour quarks to form spin-0 
pairs ("diquarks"). For a review of such effects, see [23]. Would such effects 
be important, the QGP should rather be treated as a boson gas, or a boson- 
fermion mixture. Such non-perturbative phenomena cannot be rigorously 
treated within QCD. 
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Figure 5.1. A typical result from eq. (5.13), showing the relation between the 
critical values of the temperature T and the chemical potential \i. The bag 
pressure B 1 ^ = 145 MeV was used. 

5.2 The formation of a QGP - experiments 

There are only two known situations where a QGP object could form 

• In a very hot and dense astrophysical environment, for instance, in the 
early Universe or in the interior of a heavy neutron star or a collapsing 
star, as will be discussed in a separate chapter 

• In a high-energy heavy-ion collision taking place in an accelerator or 
due to cosmic rays. 

A central high-energy collision between heavy nuclei, resulting in a compres- 
sion of nuclear matter, is the most promising process fir creating a QGP 
in a laboratory environment. There are presently several such experiments 
running at the SPS collider at CERN and the AGS collider at Brookhaven. 
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Two more colliders are under construction, LHC at CERN and RHIC at 
Brookhaven. 



Facility 


Location 


Starting date 


Energy / nucleon 


AGS 


Brookhaven, USA 


1986 


4.84 GeV 


SPS 


CERN, Switzerland 


1986 


17.2 GeV 


RHIC 


Brookhaven, USA 


1999 


200 GeV 


LHC 


CERN, Switzerland 


~ 2005 


5.4 TeV 



There are three distinguishable processes in models of ultra-relativistic heavy- 
ion collisions and QGP formation. The first is the collision of the two nuclei 
and the formation of the thermalized QGP. Several models have been devel- 
oped to describe this dynamical process, e.g. the partonic cascade model [24] 
and the QCD string-breaking model based on the Lund string model [25]. An 
overview is given in [26]. The second process involves the (relativistic) hy- 
drodynamic behaviour of the plasma. The third process is when the plasma 
has cooled to the critical temperature, T c m 150 — 200 MeV, and starts to 
hadronize. 

Several detailed calculations have been presented [27-29] in support of the 
view that a thermalized QGP will be produced in the planned experiments 
at RHIC and LHC. 

5.3 The decay of a QGP 

The decay of a QGP is believed to take place along one of two possible decay 
channels. It may expand until the density drops below the threshold of 
stability for quark matter, or it may hadronize through emission of particles, 
mainly pions, from the surface of the confining bag. 

It can be shown by a fairly simple reasoning [4] that pion evaporation 
cannot be the dominant process in the decay of the plasma. There are three 
possible processes that can create pions on the surface of a QGP bag: 

• A quark loses kinetic energy while bouncing off the surface from inside. 

• An antiquark loses kinetic energy while bouncing off the surface. 

• A quark and an antiquark annihilate at the surface, emitting a pion. 
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The energy emission through these three processes is [4] 



d 3 E T 6 , N 

0.02— (5.15) 



dtd 2 x f% 

where f n is the decay constant of the pion and x is a spatial dimension. 
The cooling by pion emission increases with temperature, but integrating 
eq. (5.15), the pion radiation falls below the thermal emission rate up to 
temperatures of 300 MeV. This leads to the conclusion that cooling by pion 
radiation is probably not a major channel of energy loss. 

This supports the other alternative, i.e., that the plasma will expand until 
it reaches the critical temperature T c and then convert into a hadronic gas, 
while maintaining thermal and chemical equilibrium. A similar approach 
[30] is that the partonic reactions inside the bag change the parton density 
locally so that a hadron can be formed. Hence, the reactions split the bag 
into smaller regions, where the hadrons form due to local density fluctuations. 

The decay due to expansion and cooling hence seems to be the dominant 
mechanism, and the lifetime of a QGP should be a few times 10~ 23 s. 
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Chapter 6 

Signatures of a QGP 



In order to detect a QGP one needs some clear experimental signatures of its 
formation and/or decay. Many such signatures have been suggested in the 
literature, and some of them are reviewed here. Searching for a QGP forma- 
tion is an event-by-event search where the events detected by the detector 
can be connected to some signature of a QGP. If the critical energy density 
is reached, several types of reactions are possible, not only QGP formation. 
In experiments one tries to select events with central collisions, i.e., with a 
small impact parameter, where high mass densities are most likely to occur. 
This can be done by focusing on events with, e.g. high multiplicities, many 
nuclear fragments or protons or a high symmetry around the beam axis. 

6.1 Kinematic probes of the equation of state 

Such probes are based on the energy density, pressure and entropy of su- 
perdense hadronic matter as a function of the temperature and chemical 
potential. One looks for a rapid rise in the effective degrees of freedom, as 
expressed by the ratios ^ or ^ over a small temperature range (s being 
the entropy). If a first-order phase transition takes place, these quantities 
exhibit a discontinuity or a step-like rise [31]. 

It is impossible to directly measure temperature and energy density. In- 
stead, one measures the average transverse momentum, < pt >, and the 
rapidity 1 distribution of the produced hadrons, These quantities can be 

1 Rapidity, y, is a commonly used generalization of the velocity 
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theoretically related to the variables T, s and E. If a rapid change occurs in 
the effective number of degrees of freedom, which would happen if a QGP 
was formed, one would see an s-shaped curve in a plot of < pt > versus 
This plot serves as a phenomeno logical phase diagram, and expresses 
the simple fact that a QGP decay is expected to give many hadrons, and 
that they would get high p^ values due to the high temperature. 

<p T > 




Figure 6.1. A schematic plot of the expected behaviour of the mean transverse 
momentum versus the density of produced particles in rapidity space. 

There are other methods built on similar ideas, reviewed in [26]. 



y = arctan/ifvii) = arctan/i( — ) = - In I ) , 

pa 2 VPo — P[| / 

where the momentum four-vector, p M = (po,P\\,P±) and po = E p /c. E p is the energy of the 
incoming particle. The parallel direction " ||" along the beam is taken as the spatial z axis, 
and "_L" is orthogonal to the beam. The rapidity is, unlike the velocity, a relativistically 
additive variable for a change of inertial system along the beam direction. 
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6.2 Electromagnetic probes 



The clearest signal of the QGP is probably an excess of produced lepton paris 
and photons. They probe the earliest and hottest phase of the evolution 
of the QGP fireball without being affected by final-state interactions. The 
drawback is that these signals are supposed to be rather weak and cluttered 
with signals from hadronic processes. 

Lepton pairs have for long been considered clear probes of the QGP. 
Thermal di-leptons are produced when a quark and an antiquark annihilate 
in a QGP. However, this so-called Drell-Yan lepton-pair production occurs 
also in normal hadron collisions. Recent calculations [32-35] show that lepton 
pairs from the equilibrating QGP may dominate over Drell-Yan production 
where invariant lepton masses are in the range 1 — 10 GeV. If lepton pairs can 
be measured above the expected Drell-Yan background from nucleon-nucleon 
collisions up to several GeV of invariant mass, the early thermal evolution of 
the QGP phase can be traced in a rather model independent way [36]. 

Also photons from the QGP compete with the background of several 
hadronic reactions. The hadronic radiation spectrum is not concentrated in 
a single narrow resonance. In fact, a hadron gas and a QGP in the vicinity of 
the transition temperature T c emit photon spectra of roughly equal intensity 
and similar spectral shape [37]. If the QGP temperature lies clearly above 
T c , a signal of photons from the QGP phase would be visible [36,38,39]. 

6.3 Strangeness enhancement, related probes 

The most commonly proposed signatures of restored chiral symmetry in dense 
baryon-rich hadronic matter are enhancements in strangeness and antibaryon 
production. The production of hadrons containing strange quarks is normally 
suppressed in nuclear reactions [40], because s quarks do not exist in the 
original heavy ions. When a QGP is formed, the production of hadrons 
containing s quarks is expected to be saturated thanks to the production 
of ss pairs [41]. The yield of strange baryons and antibaryons is therefore 
predicted to be strongly enhanced in the presence of a QGP [42,43]. 

There has been many observations of elevated strange-particle production 
in nuclear collisions. However, such an enhancement alone cannot prove the 
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presence of a QGP. Strange particles can be copiously produced in hadronic 
reactions, and several calculations have been presented [44,45]. Nevertheless, 
the enhancement of A production over a wide rapidity range, observed in 
heavy-ion collisions at 200 GeV/nucleon does not seem to be explained by 
these models. Unfortunately, it would be premature to conclude that a QGP 
has been formed. The QGP calculations use the assumption that final-state 
interactions in the decay of the QGP do not modify the hadron yield, while 
there are strong reasons to believe that the strangeness-carrying hadrons 
interact strongly. Therefore final-state interactions could very well erase all 
traces of the plasma phase. 

A strangeness abundance can still be a useful trigger for verifying other 
signals of a QGP. 

As stated earlier, a phase transition from hadronic matter to QGP would 
restore the broken chiral symmetry, i.e., make the quarks behave as if mass- 
less, or very light. This symmetry breaks down again when the plasma decays 
into normal hadrons resulting in the formation of so-called disoriented chi- 
ral condensates, DCC [46]. The DCC would later decay into neutral and 
charged pions, favouring a neutral-to-charge ratio different from the value of 
1/2 from isospin symmetry. Such events that virtually violate isospin sym- 
metry have been observed, namely the so-called Centauro events from cosmic 
rays interacting in the atmosphere [47]. 

6.4 Hard QCD probes of deconfinement 

It is commonly believed that a cc pair produced by fusion of two gluons 
cannot bind inside the QGP. Therefore, J/ip production is suppressed in 
the collision of two nuclei forming a QGP compared to a normal hadronic 
process. This assumption is based on the fact that a bound state cannot 
exist when the colour screening length Ad ~ ^ is less than the bound- 
state radius < > x / 2 [48]. Computer simulations have shown that this 
seems to happen slightly above the transition temperature T c , say at T ~ 1.2 
T c . The formation of the cc system takes a time of about 1 fm/c, and 
therefore it could require a rather long QGP lifetime before the suppression 
becomes effective. The details of modeling J/ip suppression near T c are quite 
complicated and the phenomenon has been extensively studied by several 
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authors [49-57,45,44]. There has been experimental observations of J/ip 
suppression in Pb-Pb collisions by the NA38 experiment at the CERN-SPS 
[58]. There are, however, other hadronic mechanisms that could explain the 
results equally well [59]. Quarkonium suppression is nevertheless believed to 
be the most promising signal so far of the formation of a QGP. 

Hard partons, i.e., quarks or gluons with very high energies, are believed 
to be formed at the early phase transition. If such a parton travels through 
a deconfined medium, it finds much harder gluons to interact with than it 
would in a confined medium where the gluons are constrained by the hadronic 
parton distribution [60]. The mechanisms in this case are similar to those 
responsible for the electromagnetic energy loss of a fast charged particle in 
matter. Energy may be lost either by excitation of the penetrated media 
or by radiation. Fast partons lose much more energy per unit length in a 
QGP than in hadronic matter. Hence the energy loss could tell whether 
or not the parton has travelled through a QGP. On the other hand, the 
average transverse momentum, < >, of produced hadrons is supposed to 
be enhanced from a QGP, due to the high temperature. A suppression of 
hard partons can therefore be difficult to disentangle from the rise in < p T >. 

6.5 Signatures - a summary 

It is obvious that a QGP can only be detected by the use of a combination 
of signals from different stages of the high-temperature phase of QCD. The 
signals discussed here all have normal hadronic counterparts. However, the 
overall signals of a QGP would presumably not be simultaneously duplicated 
by normal hadronic reactions. 

The QGP has yet to be uniquely observed or identified, but there are 
data that appear as rather promising. One such example is the detection 
of an enhancement of strange particles made at both the AGS [61] and the 
SPS [62-66] . The enhancement in the A yield measured over a large rapidity 
interval is difficult to describe from purely hadronic interactions. An equi- 
librated QGP with strangeness neutrality and large strangeness saturation, 
which hadronizes and decays instantaneously, fits the SPS data [67-69], but 
further experimental information is necessary to rule out other explanations. 
The current (1996), experimental situation is discussed in [70]. 
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Chapter 7 

Astrophysical aspects of QGP 



There are two situations in astrophysics and cosmology where the QGP phase 
could be relevant. The first is inside compact stars, i.e., neutron stars and, 
possibly, so-called quark stars. The second situation is in the early Uni- 
verse, before the quark-hadron transition. There are other areas connected 
to those, which might involve a QGP (one interesting example being during 
a supernova collapse) but these are of a more speculative nature and will not 
be dealt with here. 



7.1 QGP in the early Universe 

The early stages of the evolution of the Universe are described by the stan- 
dard Big Bang scenario. Its empirical basis are the following observations: 

• The redshift in the spectra of distant galaxies due to the expansion of 
the Universe. 

• The existence of the 2.7 K cosmic background radiation. 

• The abundance of light elements, especially the He/H ratio. 

• An anisotropy in the background radiation corresponding to a temper- 
ature fluctuation, AT/T « 6 • KT 6 . 

Indirect information about the early Universe is in principle obtainable either 
through electromagnetic or gravitational radiation, or, on the more specula- 
tive side, through exotic relics like very massive particles or small black holes. 
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In practice, the only observations so far are from electromagnetic radiation. 
This radiation gives information about the Universe when it was older than 
~ 2 x 10 5 years. This is due to the fact that at that time the temperature of 
the Universe was somewhat lower than the corresponding binding energies 
of electrons in light atoms. The electrons and the light nuclei then formed 
stable atoms and the Universe became transparent to photons. The cosmic 
background radiation originates from this period. 

The abundance of light elements supports the standard model of primor- 
dial nucleosynthesis, i.e., the formation of light nuclei directly after the Big 
Bang. This process took place when the Universe was a few minutes old, and 
this is the earliest epoch of which there is more or less certain information. 

If one wants to study the very young Universe, i.e., when it was younger 
than a few minutes, an extrapolation of the cosmological model to earlier 
times is required. The basic physical theory used in cosmology is Einstein's 
general relativity, but one also needs a model for high-energy particle inter- 
actions. 

Hadrons were formed in the quark-hadron transition when the Universe 
was about 10 -5 s old. The temperature was then about 200 MeV. Before 
this transition, all matter in the Universe was contained in a small region 
with an enormous density, a QGP. The transition is assumed to have been 
of first-order but this has not been shown within QCD. 

Further back in time, about 10~ n s after the Big Bang, the electroweak 
phase transition took place. The temperature was then about 100 GeV, 
i.e., in the energy range of modern particle accelerators. At T ~ 100 GeV 
there is a spontaneous symmetry breaking 1 of the electroweak SU(2) xU(l) 
symmetry, through a weak first-order phase transition. All leptons then 
acquire mass, while the intermediate vector bosons split up into the massive 
W ± and Z° and the massless photon. 

Even earlier, at t ~ 10~ 37 s, when T ~ 10 15 GeV, another spontaneous 
symmetry breaking took place. Here, a unified strong and electroweak inter- 
action is believed to have been split up by the Higgs boson. 

^ee Appendix B for the role of spontaneous symmetry breaking in phase transitions. 



39 





-40 


-30 


-20 


10 





10 




20 


Time 




















► 

log(fs) 




18 


13 


8 3 




-2 


-7 


-1 


2 


Temp. 




















► 

log(T/GeV) 


Planck scale, 




Electroweak 


Nucleo- 


Now 




quantum gravity 




phase transition 


synthesis 









GUT phase 
transition 



Quark-hadron Cosmic 
phase transition background 
radiation 



Figure 7.1. Thermal history of the Universe. 

At t ~ 10 -43 s when the temperature was T ~ 10 19 GeV, (the "Planck" 
temperature), quantum gravity effects were still important. At this time, the 
cosmic horizon was about the size of one particle, and since all mass/energy 
in the Universe was contained in this small volume, gravity was important. 
Hence, a theory of quantum gravity is needed to accurately describe this 
period. All efforts to construct such a model have, however, failed. 

The fact that a quark-hadron transition took place in the early Universe 
is quite undisputed, although one may dwell about whether or not the quark 
phase was a QGP according to the standard particle physics description. The 
main reason for the recent interest in this transition is the implication the 
exact nature of the transition could have on the abundance of light elements 
in the Universe [71]. However, it has also been argued [72] that the exact 
QGP behaviour is not so crucial for the nucleosynthesis. 

The actual transition is believed to have taken place by hadronic bubble 
growth inside the QGP. This was suggested by Witten [73], and although 
some of his conclusions have been criticized [74,75], the bubble nucleation 
mechanism is still in favour. 

Witten suggested that at a temperature just below T c , bubbles of hadronic 
matter appear. Since this corresponds to a first-order transition, a difference 
in energy density between the different phases appeared; a so-called latent 
heat. As such bubbles expand they expel heat, and a pressure equilibrium 
between the two phases makes it possible for them to coexist. 

When the Universe expanded the hadronic bubbles also expanded and 
hence the formation of hadronic matter continued. When about 50% of the 
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QGP 



Figure 7.2. Bubbles of hadronic gas (HG) begin to form inside the QGP. 

total volume had converted into hadronic matter, the QGP phase began to 
form bubbles inside a hadronic sea. The further expansion of the Universe 




Figure 7.3. Here more than 50% of the total volume is occupied by hadronic 
gas, and the QGP shrinks to bubbles. 

resulted in a loss of heat, while the QGP bubbles shrank until they finally 
disappeared. Thermal equilibrium seems to be a fair assumption for these 
rather slow processes [73]. Witten also suggested the primordial QGP might 
have survived the quark-hadron transition and constitute some of the dark 
matter in the present Universe. 

For a more thorough description of the present understanding of quark- 
hadron transition in the early Universe, see [76]. 
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Figure 7.4. The further expansion of the Universe made the QGP bubbles 
shrink until they disappeared and led only the hadronic gas to participate in 
the primordial nucleosynthesis. 

7.2 QGP and compact stars 

This section reviews the possible existence of QGP in compact stars, i.e., 
neutron stars and so-called quark stars. The material is to a large extent 
taken from Glendenning's book "Compact Stars" [77]. 

The birthplace of a star is a cloud of interstellar gas. The cloud consists 
mostly of molecular hydrogen with varying density. The temperature in the 
cloud varies from 10 K to 2000 K in different spatial regions. When a local 
cluster of high-density regions experiences a perturbation, e.g. a shock wave, 
it will collapse gravitationally if the mass of the cloud is near the critical 
" Jeans" mass. The compression continues until the internal thermal pressure 
balances the gravitational contraction energy. The core has then reached a 
high-enough temperature and density to fuse hydrogen into helium. The 
contracted cloud has become a main-sequence star. 

The thermonuclear fusion continues even after all hydrogen in the core has 
been used. Now helium burns into carbon, and the fusion process continues 
until iron has been formed. 

These processes do not depend on the mass of the star. The fusion rate 
is faster for heavy stars, but the processes are the same for all stars. 

It is important to remember that fusion will start when and where the 
temperature and density is high enough in the star. At the beginning, fusion 
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takes place only in the core, but the temperature and density in the outer 
layers of the star may be high enough for fusion to spread. 

When fusion stops in the core of a light star of a few solar masses M , 
the thermal pressure decreases and the star starts to contract again. This 
increases the temperature and density all over the star, making hydrogen 
burn in the outer layers of the star. Such a star is called a red giant. All this 
repeats for helium, carbon etc. These processes can occur in an explosive 
manner, ripping away most of the star to form a planetary nebula. The 
remaining core of the star is too light to maintain fusion and forms a white 
dwarf with a surface temperature of around 8000 K. It radiates energy for 
about 10 10 years and thereafter becomes a black dwarf. 

For a heavier star (m > 8 M Q ) the evolution is somewhat different from 
the period when iron has formed in the core. The burning process then 
continues in the outer layers of the star, adding to the iron content in the 
core. Gravity compresses the core until the electrons are captured by the 
protons; a process called inverse beta decay. The pressure in the core sud- 
denly decreases and an enormous implosion occurs, which raises the central 
temperature in the core up to several tens of MeV (~ 10 11 K). The core 
neutronizes, i.e., turns electrons and protons into neutrons. It decreases in 
volume, creating a shockwave that rips the outer regions of the star apart 
and creates a supernova in the sky. The total energy release is some 10 46 J. 
At this point the core has reached its final equilibrium state, composed of 
neutrons, protons, hyperons and leptons, and at the very centre, maybe a 
QGP. 

If the mass of the star exceeds the so-called Oppenheimer-Volkoff mass 
limit, where the internal thermal pressure cannot balance the gravitational 
compression, a black hole is formed. There are other types of collapse sug- 
gested for the formation of a black hole. See [77] for a review. 

The density in the core of a neutron star is enormous (~ 10 15 g/cm 3 ), 
which means that the energy density inside a neutron star could be high 
enough for a QGP to form. Since the core has a higher density than the outer 
regions, a QGP core should be surrounded by a hadronic crust. Such stars 
are called hybrid stars. In the extreme case of a so-called quark star there is 
no hadron crust at all. For further reference, see chapter 8 in Glendenning's 
book [77]. 
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Two possibilities exist concerning the existence of a QGP inside a hybrid 
star. It could have a pure QGP core with a distinct phase boundary against 
the hadronic crust. But there could also be a mixed phase with QGP and 
hadronic matter. Such a mixture could even be in the form of a crystalline 
lattice of various geometries, with the rarer phase immersed in the dominant 
one. Fig. 7.5 shows the internal structure of a hybrid star. 




Figure 7.5. A section through a neutron star that contains an inner region of 
mixed phases including a crystalline lattice of deconfined and confined matter 
and a surrounding liquid of neutrons and leptons. A thin crust of heavy ions 
forms the stellar surface [77]. 

Is it possible to model the internal structure of a hybrid star? With a 
bag constant B l l A = 180 MeV, and the coupling constants lying in the range 
of accepted values, the population of quarks, baryons and leptons will be as 
in Fig. 7.6. 

A neutron star is normally too cold to be directly observed, with one 
exception though; so-called pulsars. A pulsar is a highly magnetized neutron 
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Figure 7.6. Population densities per unit volume in a hybrid star with mass 
I.42 M® versus distance from the star centre [77]. 

star that rotates very quickly and thus emits radiation along its axis of ro- 
tation, creating a pulsed signal with a very stable period. This was actually 
discovered in 1968 by Bell [78], and the total number of discovered pulsars 
is now around 700. They are observed under various circumstances, usually 
as isolated sources, but sometimes in binary orbit with another white dwarf 
or neutron star. 

Pulsars lie in the ms to s range with an average period of 0.7 s. They are 
very stable and the fastest pulsar known, the PSR 1913+16, has a period 
of 1.55780644887275 ± 0.00000000000003 ms. The accuracy here is due to 
terrestrial time standards. The pulsars are therefore the most accurate clocks 
in the Universe. 

As one might guess, there is a physical limit of the rotational velocity 
for a neutron star. When the centrifugal forces overcome the gravitational 
contraction, the star will break up and emit matter. For a neutron star of 
mass 1.44 M this limit corresponds to a period of 0.3 ms (assuming that the 
structure of the star is optimized to achieve maximum rotational frequency). 

As can be seen in Fig. 7.7, the central density required to achieve submil- 
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Figure 7.7. The least possible theoretical central density of a star with mass 
1.44 M & required for stable rotation versus the rotational period. The energy 
density e c is normalized to e^, the normal nuclear energy density [77]. 

lisecond periods grows rapidly with decreasing period. One can thus conclude 
that matter under normal nuclear density is not able to keep the star together 
if the rotational period is of the order of 0.1 ms. A reasonable assumption is 
therefore that if a submillisecond pulsar is discovered, its core must be denser 
than neutron matter, i.e., containing QGP in some form. 

7.3 QGP and dark matter 

The concept of cosmic dark matter is well established today. There are 
clear signals from different observational sources for discrepancies between 
the observed mass density and the dynamical behaviour of galaxies. The 
luminous mass simply cannot account for the observed rotations. This has 
been known for a long time, and many models have been proposed in terms 
of non-luminous "dark" matter, which ultimately could account for 90% of 
the total mass of the Universe. 
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One such model was proposed by Witten [73] in 1984. The dark matter 
was suggested to be leftover quark objects from the quark-hadron transition 
in the early Universe. Witten showed that under certain circumstances, QGP 
bags, or " quark nuggets" as he called them, could survive to the present and 
form lumps of quark matter of size ~ 10 cm with a density of ~ 10 15 g/cm 3 . 

There are more recent, and detailed, analyses of Witten's idea. Madsen 
argued [79] that QGP bags with baryon number A< 10 35 cannot be dark- 
matter candidates. Alan, Raha and Sinha [80] claimed that QGP bags with 
baryon number 10 39 < A < 10 49 are cosmo logically stable and could solve 
the dark matter problem and even close the Universe gravitationally. 

Explaining dark matter with QGP objects is appealing in the sense that 
it relies on the well-established idea that all nuclear matter was once in the 
form of QGP. 
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Chapter 8 

Stability of cosmic QGP objects 



8.1 Introduction 

The second part of this thesis is an investigation of the possibility of stable 
QGP bags. These bags are assumed to be leftovers from the quark-hadron 
transition in the early Universe, following Witten's idea [73]. The basic idea 
is that if a pressure equilibrium occurs at the edge of the QGP bag, the 
bags could be stable. The dominant decay mechanism for a QGP is cooling 
by expansion and then hadronization, but if a pressure equilibrium would 
be present at the edge of the bag, no expansion would take place and the 
QGP would be stable. This possibility opens up new ways of explaining 
astrophysical phenomena such as dark matter and gamma-ray bursts. 

This investigation can be divided into two parts: 

• To determine the size where the gravitational pressure balances the 
other pressures occurring at the edge of the bag. 

• To determine the corresponding QGP mass and compare the result 
with the Schwarzschild radius. 

The calculations necessary to achieve these goals will be performed in several 
different configurations, varying the equation of state in different ways. 

The question of quark matter stable against gravitational collapse have 
been investigated before, in the sense of stable quark or hybrid stars, e.g. 



48 



[84,85]. Compared to those models, the present work is built on similar 
equations of state, but with a different formation mechanism, since it assumes 
that the QGP objects are primordial. 

8.2 The model 

The model used for the stable QGP bag is simple. Four different pressures 
exist: 

• A gravitational pressure. 

• A vacuum pressure, characterized by the bag constant B. 

• A kinetic pressure, from the motion of the quarks. 

• A Fermi pressure of the quarks. 

The kinetic pressure is identically zero if the temperature is zero. If a 
stable solution is to be found some kind of repulsive force must be taken into 
account. This is the degeneracy pressure of the quarks. The total pressure 
in the QGP bag will be divided into three parts: 

Pt=Pdg+Pk- B. (8.1) 

B is the negative (inward acting) vacuum pressure and pk is the outward 
kinetic pressure. pd g stands for the sum of the gravitational and the degen- 
erate (Fermi) pressure. The total pressure p t will start at some value p c at 
the centre of the bag and decline to zero at the edge of the bag. The actual 
calculations rest upon certain approximations and assumptions: 

• The quarks in the QGP bag are u, d and s quarks, all taken to be 
massless. 

• The plasma is taken to be a perfect fluid, i.e., the pressure is isotropic 
in the rest frame of each fluid element so that shear stress and heat 
transport can be neglected. 

• The plasma is assumed to behave like a mixture of relativistic boson 
and fermion gases (gluons and quarks). 
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Figure 8.1. Phenomenological picture of the pressures acting on the QGP 
bag. 

• The quark-hadron transition in the early Universe was a first-order 
transition. 

• The metric used to describe a static spherical star can also be used to 
describe a QGP bag, i.e., the so-called Tolman-Oppenheimer-Volkoff 
eq. (9.1) is valid for a QGP bag. 

• One of the parameters needed to be set at each calculation is the cen- 
tral pressure, p c . The interesting range of this parameter is assumed 
to begin at the pressure equivalent to normal nuclear energy density, 
£nuc ~ 170 MeV/fm 3 (corresponding to ~ 3 • 10 14 g/cm 3 ). 

• No decay process other than cooling by expansion and later hadroniza- 
tion is taken into account. 
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Chapter 9 



The crucial equations 

9.1 The Tolman-Oppenheimer-Volkoff equa- 



The standard equation describing the interior of a compact star is the Tolman- 
Oppenheimer-Volkoff (TOV) equation [86]. It is derived from the Einstein 
field equations with a static, spherically symmetric metric governing the grav- 
itational field for a static spherical star. In addition, the mass-energy tensor 
for the system is taken to represent a static perfect fluid. In the original cal- 
culations, the energy density and the pressure were taken to be zero outside 
the star. This leads to a boundary matching where the mass, integrated from 
the centre of the star to the edge, is set to be equal to the mass term in the 
external Schwarzschild metric. 

However, when applying this formalism to an MIT bag, with its external 
vacuum pressure B, the energy density outside the bag will be 4B. This 
seemingly strange assumption, inherent to the MIT bag model, does not 
significantly influence the metric. 

The TOV equation is derived from the assumption of a static spherical 
star taken as a perfect fluid. It is derived in many textbooks in general 
relativity, e.g. in [87]. It results in the following set of equations: 



tion 



dp 
dr 



[e(r) + p(r)][m(r) + 4nr 3 p(r)} 
r[r — 2m(r)] 



(9.1) 
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^1 = 4vrr 2 e(r), m(r) = 4tt f e(r')r' 2 dr' ', (9.2) 
dr Jo 

p(r = 0) = p c , (9.3) 
p( r = R) = 0. (9.4) 

These equations are expressed in gravitational and natural units, c = h = 
G — 1. p is the total pressure, e is the energy density, m is the mass and r 
is the radial coordinate. 

The TOV equation requires just a static spherical mass-energy distribu- 
tion and a perfect fluid. It does not depend on a certain equation of state 
that relates energy density and pressure. 



9.2 The equation of state 

The equation of state for a non-interacting QGP consisting solely by u and 
d quarks along with gluons was discussed in section 5.1. Allowing also s 
quarks changes the equation of state somewhat. The pressure is, as previously 
stated, divided into four parts, where the kinetic part is derived for a mix 
between three fermion gases and one boson gas. These circumstances change 
the equation of state to 

e(r) = 3p t (r) + AB. (9.5) 
The kinetic pressure is now with three quark flavours (/): 

and the chemical potential /j,f is given for each quark flavour. When taking 
the lowest-order gluon interaction into account, the kinetic pressure changes 
into 

87T 2 . ( 15o; s \ T 7 9 . / 50a: s \ . 

Pk = T 1 + > — vr 2 T 4 1 + 9.7 

y 45 V 47r / f L 60 V 2l7r / 

This expression has been derived in [88-90]. Here one can obviously study 
also the case when quarks do not interact, by setting a s = 0. 
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After introducing p t into the TOV equation one gets 



dpt_ _ (4p f + 4g) [4tt / r (3gt + 4g) r /2 dr ; + 4%r 3 Pt ] 

dr ~ r [r - 2 • 4tt /J" (3p t + 45) r /2 dr'] ^ ' ' 

This integro-differential equation cannot be solved analytically, and basic 
numerical methods will be used in the following calculations. 



Chapter 10 
The method 



10.1 The algorithm 

Equation (9.8) is solved with two simple numerical algorithms. The deriva- 
tive is changed according to 



dp Pi+i - Pi 



;io.i) 



dr h 
and the integral is changed into 

4tt / (3p t + 4B)r' 2 dr' ^ ■ {[p t (0) + 4B] ■ 2 + 2[p t (h) + (10.2) 
Jo ^ 

+AB] ■ h 2 + 2\p t (2h) + AB] ■ (2h) 2 + ... + \p t (nh = r) + AB] ■ r 2 } , 

where h is the constant step length. rThese discretizations change the TOV 
equation into 



P^+l ~ Pi _ (4j>i + 4g) (Emt +4^) 

n(^-2E mt ) 

where the integral approximation is denoted X)mt- 



10.3) 



As can be seen, there is a singularity at r« = 0. This makes it impossible 
to calculate p(r = h) from eq. (9.8). This problem is not new, and methods 
have been developed to cope with it. A polynomial expansion will be used 
for the first step, i.e., in calculating p(r = h), so that p(r) = YljPj rJ ' ■ Hence 
near r = 0, 

P = P{tc) + (p c T c /e c ) ( e - e c ) + ... . (10.4) 
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T c is the "adiabatic index" d[ln(p)]/d[ln(e)] evaluated at e = e c , which also 
can be put on the form dp/de taken at the same energy density e c . This 
expansion can now be written as 

p(r) ^p c -2n (e c + p c ) (p c + ^e^j r 2 + (10.5) 

so that p(0) = pc, and p(h) can be calculated. The resulting formula 

M4 P , + 4g) (£„„ M) fl06l 

P,+1_P - r,(r,-2E„J (1 °' 6) 

is to be iterated until the pressure p t = is reached at some radius r, con- 
taining the full mass (%2 int )- 

A scaling procedure has been used. All involved quantities are in km. 
This is possible by using gravitational and natural units: 



[p] 


= 1/km 2 , 


[e] 


= 1/km 2 , 


[T] 


= 1/Vkm", 




= 1/Vkm, 


[r] 


= km, 


[B] 


= 1/km 2 , 


[ml 


= km. 



The following relationships hold: 

B 1 / 4 , T, N : 1 MeV = 1 • ^ ' ) = 6.4421 • 10~ 4 1/vS, (10.7) 

2 fill 5 • 10 -4 

e,p : 1 MeV/fm 3 = 1 • = 1.3234 • 10~ 6 1/km 2 , (10.8) 

197.33 

1 4766 

m : 1 kg = 1 ■ 1,989-1030 = 7 - 4238 ' 10_31 km - ( 10 - 9 ) 
When using the scaling expressed above, the normal nuclear density e nuc = 
3 • 10 14 g/cm 3 « 170 MeV/fm 3 becomes e nuc « 2.2 • 10~ 4 1/km 2 and bag 
constant B = (150) 4 (MeV) 4 becomes B « 8, 7 • 10~ 5 1/km 2 . 
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10.2 Numerical stability and errors 



The difference formula (10.1) for the derivative has an error of order 0(h), 
while that of the integral approximation has an error of order 0(h 2 ). It is in 
principle possible to reduce the step size h, but it is then possible that one 
has to use the series solution, eq. (10.5), to determine the pressure close to 
the singularity, i.e., more than one step out of r = 0. 

When it comes to the stability criterion of the iteration procedure no the- 
oretical examination has been made. The approach here has been to reduce 
the step size by 50% as long as significant differences, (> 1%), remained 
between two iterations in the mass/radius relationship for the QGP bag. 
The final step size was 10 m, i.e., h = 0.01, which proved to be a suitable 
compromise between accuracy and CPU time. 

10.3 Parameter variations 

In the calculations there are seven variable parameters; B, a s , T, /i u , /i d , /i s 
and p c . The central pressure has been varied between 1.7 — 1.3 • 10 5 times 
the corresponding normal nuclear density. The temperature has been varied 
between and 2000 MeV, (10 12 K). No significant change in the bag size was 
caused by this variation of temperature. The same is true for the chemical 
potentials /x n , Hd and fi s . The strength of the QCD interaction, a s has been 
varied between and 2.5, which did not alter the bag size significantly. The 
only important variation that did alter the bag size (apart from that of the 
central pressure) was in the bag constant. The value of was varied 
between 120 and 180 MeV. 
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Chapter 11 

Results and discussion 



11.1 Results 

When calculating the stable configurations for these QGP objects, a stability 
criterion built into the TOV equation must be considered. This is known as 
Le Chatelier's principle and can be formulated as [77] 

dM(e c ) 



9e c >0 - (1L1) 

When applying this condition to the results presented here, the number of 
stable configurations is limited. As shown in Fig. 11.1, three regions in the 
range of the applied central energy density e c allow stable configurations. 
The first region is only interesting in the upper part, since at lower e c , a 
QGP will not form. In the second region, e c is far above the QGP transition 
allowing stable configurations to occur. One should notice that at these e c 
values, the chemical potential for the c quark is the same order of magnitude 
as its mass and therefore, the c quark could be present in the centre of the 
QGP. This implication is not included in the calculations. The small third 
region, represents central densities about a thousand times higher than those 
in the Universe at the time of the quark-hadron transition. 

Presenting the integration of the TOV equation in a mass-radius diagram 
normally results in a spiral-like curve. This is true also for the formalism 
presented here, as shown in Fig. 11.2. 

If a QGP bag is to be stable, i.e., to be able to exist under a cosmological 
period of time, the temperature should be zero so that no thermal photons 
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Figure 11.1. Mass vs. applied central energy density. Stable configurations 
occur only in the density intervals with the full-drawn line. The parameter 
values are B 1 / 4 = 150 MeV, T = MeV, fid, /x s = and a s = 0. 
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Figure 11.2. This plot shows stable (full line), and unstable (hatched line) 
configurations of the QGP bag. The parameter values are B 1 ^ = 150 MeV, 
T = 200 MeV, fi u , /i^, \l s — and a s = 0. The central pressure is varied from 
10~ 5 to 100 km~ 2 , corresponding to a central energy density of 1.7 to 1.3- 10 6 
times s nuc . 
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Figure 11.3. The mass-radius relation for varying bag pressure B (full line): 
B l l A =180 MeV (curve 1.), 150 MeV (2.), 120 MeV (3.), 100 MeV (I), 
75 MeV (5.) and 50 MeV (6.). The hatched line corresponds to the 
Schwarzschild limit. All other parameter values are the same as in Fig. 11.2. 

are emitted. If the temperature is set to zero in the calculations no visible 
change in the mass-radius relationship occurs within the accuracy of the 
iterative procedure, i.e., with a step in r of 10 m. 

The eq. (9.7) for the kinetic pressure pk does not depend much on the 
strong coupling constant a s , even for a s values up to 2. Consequently, the 
mass-radius relationship is, in practice, independent of a s . The same result 
is achieved when varying the chemical potentials of the quarks. The only 
parameter variation that influences the mass-radius relationship substantially 
is the bag constant as can be seen in Fig. 11.3. 

With the choice of parameter values in Fig. 11.2, it is obvious that one 
has two possible maxima regarding the mass of the bag. One lies around 10.5 
km with a mass about 1.9 M Q . This is achieved when the applied central 
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Figure 11.4. The same results as in Fig. 11.2 complemented with the mass- 
radius relation corresponding to the Schwarzs child radius (upper hatched 
line). 

energy density is about 10 e nuc . The other lies at 8.5 km with a mass about 
1.4 M . Here, the applied central energy density is roughly 10 4 e nuc . 

Using the basic formula for the Schwarzschild mass-radius relationship, 
r s = 2Gm/c 2 , it is improbable with a gravitational collapse into a black hole, 
even at maximum mass configurations. One might ask if there is any internal 
part that has a mass-radius relationship near that given by the Schwarzschild 
relation. One then needs to look at the energy density profile of the bag and 
compare the integrated mass with the Schwarzschild mass. No such situation 
appears with reasonable parameter values. 

When varying the bag pressure B, as done in Fig. 11.3, the regions in e c 
resulting in stable mass-radius configurations is slightly shifted. One can also 
see, that significant changes in B affect the mass-radius configurations, but 
not in a way that makes gravitational collapse more probable. The value of 
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B is traditionally taken from quantum-mechanical calculations involving the 
MIT bag model, where B has been adjusted to make the calculated hadron 
masses agree with experimental data. In this thesis, the same values are 
used in a vast extrapolation to macroscopic objects, which is customary in 
astrophysical applications of the MIT bag model. However, one cannot be 
certain that the value of B is the same regardless of external circumstances. 
In the calculations presented here, it is assumed that B does not alter with 
the outside metric. The possibility that non-Euclidian geometry change the 
value of B is not examined. 

11.2 Discussion 

Based on these calculations, one can conclude that a QGP bag obeying the 
chosen equation of state has a radius of either 7.5 - 8.5 km with mass 1.2 - 
1.4 M or 9 - 10 km with mass 1.6 - 1.8 M Q , while still being stable against 
evaporation and gravitational collapse. This result agrees with other quark 
matter stability calculations [84,91]. The size and mass of these bags make 
them interesting as dark matter candidates as well as possible sources for 
gamma-ray bursts [92]. 

The calculations also show that the connection between e c and the average 
energy density is rather weak, at least when e c > 100e nuc . At high e c , the 
energy density drops very fast with increasing radius close to the centre and 
assumes a behaviour rather similar to other, lower e c configurations, further 
out. This fact could be due to numerical errors close to the centre and 
remains to be investigated. 

One circumstance that seriously could affect the possibility for QGP bags 
to survive the quark-hadron transition is the cosmic horizon size at the time 
of the transition, i.e., 10~ 5 s after Big Bang. At this time, [93] states that 
the radius of the Universe was about 1CT 11 times the present size, i.e., about 
10 12 km. 

The standard model in cosmology states that the average density in the 
Universe at the transition was about 10 20 kg/m 3 . It is possible that local 
density fluctuations could have created regions that contained QGP at a con- 
siderably higher density than the average at the time of the transition. These 
regions later expanded with the Universe and possibly split up into smaller 
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regions of QGP when fluctuations inside the region started local hadroniza- 
tion processes. Gravity retarded the expansion of these regions and the size 
necessary for a stable QGP bag could have been achieved. This process, 
when gravity counteracts the spatial expansion is a dynamical process that 
would be very interesting to look further into. Also, it would be of crucial 
importance to know if a QGP can exist in the intermediate density region 
between e nuc and a few times e nuc . This region is considered to lie below the 
density needed to create a QGP out of compressed nuclei. However, such 
investigations lie beyond the scope of the present thesis. 
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Appendix A 
Chiral symmetry 



This Appendix on chirality and chiral symmetry follows roughly Donoghue, 
Golowich and Holstein [81]. 

If one lets ifi(x) be a solution to the Dirac equation for a massless particle: 

% fill) = (A.l) 

and multiply this equation from the left by 75, one can use the anticommu- 
tativity with 7^ to obtain 

i M = 0. (A.2) 
The following combinations can be formed: 

^ = ^(1 + 75)^, ^=^(1-75)^ (A.3) 

The quantities ipL and ip R are solutions of fixed chirality, in this case also 
called handedness. For a massive particle moving with precise momentum, 
these solutions correspond to the spin of the particle being, respectively, an- 
tialigned (left-handed) and aligned (right-handed) relative to the momentum. 

The matrices Tlr = (l±7 5 )/2 are chirality projection operators, fulfilling 



T L + T R -1 

r L r L = r L , r R r R = r R , r L r R = r R r L = o. 1 ] 

Chirality is a natural label to use when referring to massless fermions thanks 
to its Lorentz invariance. A left-handed particle is left-handed to all ob- 
servers. A physical example of this is the neutrino in the standard model. 
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Neutrinos are left-handed chiral particles, i.e., if they are massless. There is 
no evidence for right-handed neutrinos. 

Chirality is easy to add to the Lagrangian formalism. The Lagrangian 
for a non-interacting massless fermion is 

L = kp fii) (A.5) 

or 

L = L R + L L = iip L ftil> L + iil> R fiipR. (A. 6) 

These Lagrangian densities are invariant under the global chiral phase trans- 
formation 

*Pl,r( x ) exp(ia L:R )i! L:R (x), (A. 7) 

where the phases oll,r are constant and real-valued. Hence the Lagrangian 
density for massless particles has a chiral symmetry. If we use the famous 
Noether's theorem 1 we can associate conserved particle number current den- 
sities 

JU = i>L,Rl^L,R = 0) (A.8) 

with this invariance. These chiral current densities form the basis for the 
vector current V fi (x): 

v» = r L + r R (a.9) 

and the axial vector current A^(x): 

A» = r L -r R . (a. io) 

It is now possible to construct conserved chiral charges: 

Q L , R (t) = J d 3 xJl R (x), (A.ll) 

and these represent the number operators for the chiral fields tpL,R- 

The vector charge Q = Ql + Qr is the total number operator; whereas the 
axial vector charge, Q$ = Ql — Qr, is the number operator of the difference 

1 Noether's theorem states that for any invariance of the action under a continuous 
transformation of the fields in the Lagrangian, there exists a classical charge Q, which is 
time independent and is associated with a conserved current J M ; <9 M J M = 0. 
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above. These charges simply count the sum and the difference of the left- 
handed and the right-handed particles. 

The u and the d quark have masses far below the QCD scale parameter 
A ~ 0.2 GeV which makes the assumption that m u ss m d ps quite reason- 
able. It is even a fairly good approximation to treat the s quark as massless, 
at least if the temperature (in energy units) is comparable to the mass of the 
s quark. Due to the connection between mass and chiral symmetry shown in 
eqs. (A.1-A.3), chiral symmetry is not present as long as the fermions have 
mass. 

When referring to QGP, I have previously stated that chiral symmetry 
is restored in the phase transition to QGP. This is only approximately so 
because the quarks inside the plasma are not massless. However, since the 
temperature and the scale parameter A is roughly 200 MeV, the assumption 
that m u md ~ m s ~ is fairly reasonable. The true chiral symmetry 
restoration also depends on the fact that the quark condensate expectation 
value vanishes at the phase transition. With these approximations, chiral 
symmetry is said to be restored in the QGP phase. 
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Appendix B 

The physics of phase transitions 



This description of the role of symmetry-breaking in phase transitions is 
inspired by the book "Cosmology" by Coles and Lucchin [82]. 

The transition from some disordered phase to an ordered phase is labelled 
by a lowering of the symmetry in the system, which can be characterized by 
an order parameter $. $ is zero in the most disordered phase and rises when 
more ordered phases appear. 

It is illustrative to take the solid/liquid phase transition as an example. 
The ordered phase has a low degree of (discrete) rotational symmetry, while 
the liquid gets an almost perfect rotational symmetry when the atoms move 
around at random. The order parameter <£> in this case represents the devia- 
tion of the spatial distribution of ions from the homogeneous distribution they 
have at T > T c , the melting point. If T > T c , $ = and if T < T c , $ > 0. 

When $ grows significantly, a symmetry-breaking transition has occurred. 
Such a transition can be caused either by an external influence, such as heat- 
ing, or by a change of the system itself. Transitions due to external influences 
are called "induced", while those due to internal phenomena are "sponta- 
neous". For the spontaneous symmetry-breaking processes it is convenient 
to consider the free energy of the system, F = U — TS, where U is the inter- 
nal energy, T is the temperature and S is the entropy. A condition for the 
system to be in an equilibrium state is that F is at a minimum. At T = the 
free energy coincides with the internal energy. At T > 0, whatever the form 
of U, an increase in S leads to a decrease in F, and is therefore energetically 
favoured. 
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Phase transitions can be of two kinds; first- or second-order. If $ rises 
continuously from zero when T <Tq (the transition temperature), the tran- 
sition is said to be of second-order, or continuous. If $ jumps disco ntinuously 
to a non-zero value just below Tc then the transition is of first order. Here, 
the entropy also exhibits a discontinuity at Tc- In a first-order transition, 
heat is adsorbed by the system in going from the low-temperature to the 
high-temperature phase. This heat is the latent heat, Ql = T C AS, of the 
transition, where AS is the entropy change. 



cD 




Figure B.l. Order parameter vs. temperature for a first-order (left) and a 
second-order (right) transition. 



The QGP phase transition involves a spontaneous breaking of the chiral 
symmetry, i.e., chiral symmetry is broken in hadronic matter, but restored 
in the QGP. The order of the QGP phase transition has not been uniquely 
determined. Two massless flavours would generate a second-order transition 
[20,21], while three massless flavours would give a first-order transition. In 
astrophysical calculations that involve a QGP phase transition, one often 
assumes that the transition is of first order. 
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